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^ : ABSTRACT 

o : 

■ Recent radial velocity observations have indicated that Jovian-type planets 

. can exist in moderately close binary star systems. Numerical simulations of the 

Q_|. dynamical stability of terrestrial-class planets in such environments have shown 

Q . that, in addition to their giant planets, these systems can also harbor Earth-like 

il3 . objects. In this paper, we study the late stage of terrestrial planet formation in 

. such binary-planetary systems, and present the results of the simulations of the 

>• . formation of Earth-like bodies in their habitable zones. We consider a circumpri- 

^ . mary disk of Moon- to Mars-sized objects and numerically integrate the orbits 

. of these bodies at the presence of the Jovian-type planet of the system and for 

different values of the mass, semimajor axis, and orbital eccentricity of the sec- 
ondary star. Results indicate that. Earth-like objects, with substantial amounts 
of water, can form in the habitable zone of the primary star. Simulations also 
indicate that, by transferring angular momentum from the secondary star to pro- 
toplanetary objects, the giant planet of the system plays a key role in the radial 
mixing of these bodies and the water contents of the final terrestrial planets. We 
will discuss the results of our simulation and show that the formation of habitable 
planets in binary-planetary systems is more probable in binaries with moderate 
to large perihelia. 



Subject headings: binaries: close — celestial mechanics — planetary systems — 
planets and satellites: formation — solar system: formation 
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Introduction 



For the past several years, the formation of terre strial planets in b i nary st ar systems 



has b een th e subject of res e arch by many authors. I Quint ana et al. I (|2002| ). I Quint ana 



( l2003l ). and iLissauer et al. I (120041 ) studied the interactions of planetesimals and proto- 
planetary objects around the stars of a Centauri, and showed that, terrestrial-type planets 
can form around thes e sta rs when dynamical friction is included in numerical simulations. 
Barbieri et al. I (120021 ) and iTurrini et al. I (120051 . |2006| ) also studied terrestrial planet forma- 
tion in this system and by considering gas-drag as the primary force for reducing planetes- 
imals relative velocities, showed that, it is indeed possible to form terrestrial- c lass o bjects 



around the primary of a Centauri stars. In a recent article, iQuintana et al. I (120071 ) have 



extended their simulations to wider binary systems, and identified regions of the parameter- 
space for which terrestrial planets can form around the stars of the binary. 

The studies of terrestrial planet formation in dual-star systems, as presented in the 
literature, share one common ground: The systems considered in these studies do not contain 
giant planets, and the formation of terrestrial planets has been simulated only at the presence 
of the two stars of the system. In this paper, we extend these studied to more complex 
environments and simulate the formation of terrestrial planets in binary star systems in 
which the primary star is host to a Jupiter-like planet. The purpose of our study is to 
understand how, in such binary-planetary systems, the dynamics of the stellar companion 
will affect the formation and the water contents of Earth-like objects in the habitable zon^ 
of the primary star. 

The systems of our interest are moderately close (< 40 AU) binaries. Recent detections 
of Jovian-type planets in such environments (e.g., GJ 86, cf. Els et al. 2001, 7 Cephei, 
cf. Hatzes et al. 2003) have raised questions about the formation of these objects and the 
possibility of the existence of smaller bodies in these systems. In regard to the latter, simu- 
lations of the orbital dynamics of terrestrial planets in the 7 Cephei planetary system have 
indicated that small objects can have long-term stable orbits in binary-planetary systems 
provided their orbits lie outside the influence zon^ of t he system's giant plan et and are lim- 
ited to the region between this object and its host star (jHaghighipour II2OO6I ). In this study, 
we choose systems in which the primary star has a Jovian-type body in an orbit outside its 
habitable zone. 



^The region around a star where a terrestrial-class planet can maintain liquid water on its surface. 

^The influence zone of a planetary-sized object with a mass is defined as the region between ap{\ — 
Cp) — 3Rh and Op(f -I- e^) + 3Rhi where Op is the planet's semimajor axis, Cp is its orbital eccentricity, and 
Rh ~ ap{mp/3AlY^^ represents its Hill radius. The quantity M denotes the mass of the central star. 
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Unlike the majority of the binary-planetary systems that have so far been discoverecjf], 
the relatively smaller separations of the binary components in our systems imply that the 
effect of the farther companion on the formation and dynamical evolution of planets around 
the primary star is not negligible. For instance, the perturbative effect of the secondary 
star can change the structure of the circumprimary disk and truncate it to a smaller size 
( lArtymowicz fc Lubow I Il994l ). The latter removes material that may be used in the for- 
mation of planets. For this reason, it was believed that circumstellar disks around the 
components of a close binary m ay not be large a n d massive enou g h to b e gin planet forma- 



tion. Howe yer, observatioii s bv iMathieu I (Il994j ) , lAkeson et al. I ( 119981 ) , [Rodriguez et al 



( 119981 ). and iMathieu et al. I (120001 ) have shown otherwise. These observations confirm the 
presence of potentially planet-forming environments around the components of binary stars 
and imply that planet formation around a star of a binary may be as common as around 
a single s tar. In fact, the observa tions of two well-separated disks in the binary system of 
L1551 by iRodriguez et al. I (119981 ) indicate that, despite of disk truncation, it is still pos- 
sible for the both components of a binary to retain a relatively significant amount of their 
original circumstellar materials (0.03 to 0.06 solar-masses) in disks with considerable radii 
(~ lOAU). The masses of these disks are c omparable to the min i mum solar-mas s model of 
the primordial nebula of our solar system (IWeidenschilling 1119771 : iHayashi Ill98ll ). implying 
that, planet formation in dual-star systems can begin and continue in the same fashion as 
around our Sun. 

In this paper, we base our study on the latter consideration. We assume that in a 
binary system, planetesimal formation follows similar process as around a single star, and 
giant and terrestrial planets are formed through the interactions of these objects. It is 
important to emphasize that in such systems, the stellar companion has a strong effect on the 
accretion of planetesimals and the formation of larger bodies. In general, t he perturbations 



due to the secondary star in crease the relative velocities of planetesimals (IHeppenheimer 



19781 : IWhitmire et al. I Il998l ) , which may cause their collisions to result in breakage and 



fragmentation. This object can also inhibit the formation of p rotoplanets by destabil izing 



the regions where the building blocks of these objects exist ( IWhitmire et al. I Il998l ). It 



has, however, been shown that the effect of the binary companion on increasing the relative 
velocities of planete simals may be counterbalanced by dissipative forces such as gas- drag and 
dynamical friction (IMarzari et al. Ill997l : iMarzari fc Scholl Il2000l : iThebault et al. II2006I ). As 
shown by these authors, the combined effect of gas-drag and gravitational force of the stellar 
companion results in a strong alignment of the periastra of planetesimals, which increases 
the efficiency of their accretion by reducing their relative velocities. In this paper, we assume 



See iHaghighipour 1 (|2006f ) for a complete and up-to-date list of binary-planetary systems. 
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that a Jupiter-like planet has already formed around the primary of our binary star syste 
and the interactions of planetesimals have been efficient and have resulted in the formation 
of a disk of planetary embryos (e.g., via oligarchic growth, cf. Kokubo & Ida 1998) around 
this object. 

As mentioned earlier, we focus our attention on the formation of habitable planets, 
that is. Earth-like objects in the habitable zone of the primary star. Since all life on Earth 
requires the presence of liquid water, we consider water-rich planets to be the best candidates 
for habitability, and pay close attention to the acquisition of water during the formation of 
these objects. Similar to the current models of the formation of habitable planets in our solar 
system, we assume that cometary material, if existed around the stars of a binary system, 
would provide little to no water to the terrestria l planets that might form in the habitable 



zone of the primary star. We adopt the model of iMorbidelli et al. I ( 120001 ). who argued that 
water-rich bodies originating in the solar system's asteroid belt were the p rimary source o f 
Earth's water, and simulate the late stage of terrestrial planet formation (jWetherill Ill996l ) 
by numerically integrating the orbits of a few hundred protoplanetary objects, for different 
values of the mass, semimajor axis, and orbital eccentricity of the secondary star. We assume 
an initial gradient in the water contents of protopl anets such that radial mixing is required 



to " deliver" water to \ 


2004; 


Ravmond 


2006| 



We identify the regions of the parameter-space of a binary-planetary 
system for which an Earth-like planet can form in the habitable region of the primary star. 

The outline of this paper is as follows. In § 2, we discuss the details of our model. 
Section 3 has to do with the numerical integrations of the system and the analysis of the 
results. In § 4, we study the formation of habitable planets, and in § 5, we conclude our 
study by reviewing the results and discussing their applications. 



"^The formation of gas-giant planets in a dual-star system is still a subject of research. While simulations 
by iNelson I (|200Cir ) indicate that gas-giant planet formation may not proceed through the disk instability 
mechanism around the primary of a binary star system with separation of ~ 50AU, recent simulations by 
Boss show that Jupiter-like planets can form in such environments via the gravitational ins t abilit y 



of a marginally unstable circumprimary disk. On the other hand, as shown by iThebault et al. ] ^04 



the core accretion mechanism may also be able to form giant planets around the primary of a binary star 
system. However, as the results of their simulations for planet formation in the 7 Cephei system indicate, 
the semimajor axis of the final gas-giant planet may be smaller than its observed value. 
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The Model 



As mentioned in the introduction, we would like to study the formation of terrestrial 
planets in the habitable zone of the primary of a moderately close binary-planetary system. 
We are primarily interested in understanding how the motion of the secondary star affects 
the dynamics of a disk of protoplanetary objects and the final assembly and water-contents of 
the resulted terrestrial-sized bodies. In other words, we would like to study how the process 
of habitable planet formation in a system consisting of a star, a disk of planetary embryos, 
and a giant planet will be altered if a stellar companion is introduced to the system. 

The statement above portrays a general picture of our model. To ensure the habitabil- 
ity of such a system. Earth-like objects have to form in the habitable zone of its primary 
star and maintain lon g-term stable orbits in that region. On the other hand, as shown by 
Haghighipour I (120061 ). terrestrial planets can have stable orbits only at distances close to 
the primary star and outside the influence zone of its giant planet. This requires that the 
habitable zone of the primary to be considerably closer to it than the orbit of its planetary 
companion. To satisfy this requirement, and also for the purpose of comparing habitable 
planet formation in binary-planetary systems with that around single stars, we make the 
following assumptions. 



1) We assume that the primary of our system is a Sun-like star. As indicated by lKasting. Whitmire fc Reync 
( 119931 ) . the habitable zone of such a star will extend from 0.95 AU to 1.37 AU. This is a 
conservative estimate that place s the outer boundary of the habit able zone at a distance 

where CO2 clouds start to form (jjones. Underwood, fc Sleep 1120051). The outer edge of thi s 
regi on may in fact be at larger distances. As shown by iForget fc Pierrehumbert I (119971 ). 
and iMischna et al. I (I2OOOI ). the outer boundary of the habitable zone of the Sun may be at 
approximately 2.4 AU from this star. In this study, we adopt a relatively moderate approach 
and consider a habitable zone between 0.9 AU and 1.5 AU for our primary star. 

2) We consider the giant planet of our system to be a 1 Jupiter-mass object. Since we 
would like to study how the orbital dynamics of the secondary star affects the interactions 
of planetary embryos, we assume the orbit of this planet is circular. We also consider the 
semimajor axis of this object to be at 5 AU, outside the habitable zone of the primary star. 

3) We choose the mass of the secondary star to be 0.5, 1.0, and 1.5 solar-masses. We consider 
the semimajor axis of this object to have the values of 20, 30, and 40 AU, and its eccentricity 
to be 0, 0.2, and 0.4. Since we are interested in the formation of habitable planets at the 
presence of the Jovian-type planet of the system, it is necessary to ensure that, for any 
combinati on of these paramet e rs, th e giant planet will have a long-term stable orbit . As 
shown by iHolman fc Wiegert I (119991 ). in order for a planet in a circular orbit, to be stable 
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in a binary star system, its semimajor axis cannot exceed the critical value (oc) given by 

a^/flft = (0.464 ± 0.006) + (-0.380 ± 0.010)/ib + (-0.631 ± 0.034)efe 

+ (0.586 ± Omi)fibeb + (0.150 ± 0.041)e^ + (-0.198 ± 0.074)/ibe^ . (1) 



In this equation, Of, and are the semimajor axis and orbital eccentricity of the stellar 
companion, and fif, = M2/M1, where Mi and M2 are the masses of the primary and secondary 
stars, respectively. Figure 1 shows the graph of Oc in term of the eccentricity of the binary. 
In this figure, /if, = 1 and at = 20, 30, and 40 AU. As shown here, a giant planet with a 
semimajor axis of 5 AU will not have a stable orbit in an equal-mass binary with a separation 
of 20 AU and an eccentricity of 0.2 or higher. Similar situation exists for a 30 AU binary 
with a 0.4 or larger eccentricity. We use equation (1) to identify the combinations of the 
mass and orbital parameters of the secondary star for which the giant planet of the system 
becomes unstable, and simulate the formation of habitable planets for those combinations 
of these parameters that result in a stable orbit for this object. 

4) We assume that planetary embryos have already formed in the region between the pri- 
mary and the giant planet. We consider a disk of 115 Moon-to Mars-sized bodies, with 
masses ranging from 0.01 to 0.1 Earth-masses, randomly distributed, by 3 to 6 mutual Hill 
radii, between 0.5 AU and 4 AU from the primary star. The masses of embryos increase 
with their semimaior a xes (a) and the number of their mutual Hill radii (A) as a^/^A^/^ 
(IRaymond et al. 1120041 ). The surface density of our disk model, normalized to a density of 
8.2 g/cm^ at 1 AU, varies as r , where r is the radial distance from the primary star. The 
total mass of our disk model is approximately 4 Earth-masses. Figure 2 shows the graph of 
one of such disks. 

5) We assume that the water to mass ratios o f embryos follow the current distribution of 
water in primitive asteroids of the asteroid belt (jAbe et al. 2000). That is, embryos inside 2 
AU are taken to be dry, the ones between 2 to 2.5 AU are considered to contain 1% water, and 



those b e yond 2.5 AU are assumed to have a water to mass r atio of 5% ([Raymond et al. 112004 



2005al lbl: I Raymond II2006I : [Raymond. Barnes &: Kaib II2006I ). We also consider an initial Iron 
content for each embryo. This value is obtained by interpolating between the values of the 



Iron contents of the terrestrial planets ( Lodders fc Fegley 1998 : Raymond et al. 2005al lbl) 



with a dummy value of 40% in place of Mercury becuase of its anomalously high Iron content. 



- 7- 



Numerical Simulations 



Using the N-body integration package MERCURY (IChambers Ill999l ). we numerically 
integrated the equations of motion of the planetary embryos of our disk model for different 
values of the mass, semimajor axis, and eccentricity of the secondary star. We allowed 
the protoplanetary objects to collide with one another and assumed that each collision was 
perfectly inelastic. We also assumed that no debris was generated during a collision, and 
that the effect of the energy released during an impact, on the morphology and structure of 
the colliding objects, was negligible. 

We carried out a total of 46 simulations, each with a time step of 6 day^. Figure 3 
shows the results of one of such simulations. In this figure, the separation of the binary is 30 
AU, its eccentricity is 0.2, and the mass of the secondary star (not shown in the figure) is 0.5 
solar-masses. As shown here, after 100 Myr, a terrestrial-sized object (1.1 Earth- masses), 
with substantial amount of water, is formed in the habitable zone of the primary star. The 
orbit of this object has a semimajor axis of approximately 1.18 AU and an eccentricity of 
~ 0.06. The Jupiter-sized planet of the system is shown as a big black circle. 



S imilar to the simulation s of terrestrial planet formation around single stars (IMorbidelli et al 



2000l : [Raymond et al. II2004I ). our simulations are stochastic. That is, for a given set of or- 
bital parameters of the binary companion, different initial distributions of protoplanetary 
objects produce different results. For this reason, for each set of the initial orbital parameters 
of the binary companion, we carried out simulations for three different random distributions 
of planetary embryos. Figures 4 shows the final result of such simulations for two different 
cases. The case on the left corresponds to a binary with a separation of 30 AU and a sec- 
ondary of 1.0 solar- mass in a circular orbit. The case on the right represents the results of 
simulations in a system in which the same secondary star is now in an orbit with a semimajor 
axis of 40 AU and an eccentricity of 0.2. Each simulation, from top to bottom, corresponds 
to different distribution of planetary embryos with the simulations on the same row having 
similar initial distributions of protoplanetary objects. 

In addition to the stochasticity of simulations, figure 4 also shows the relation between 
the orbital eccentricity of the stellar companion and the water contents of the final bodies. 
As shown here, for identical initial distributions of planetary embryos (i.e., simulations on 
the same rows), the total water content of the system on the left, where the secondary star 
is in a circular orbit, is higher than that of the system on the right, where the orbit of the 
secondary is eccentric. Figure 5 shows this for several other simulations. As depicted by 



^Slightly smaller than 1/20 of the orbital period of the closest embryo at 0.5 AU 
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this figure, for a given separation of the binary, the accumulative water content of the final 
planets decreases as the eccentricity of the binary becomes larger. 

The fact that the final assembly of terrestrial planets in a system with an eccentric 
secondary star contains less water implies that, prior to the formation of these objects, 
most of the water-carrying embryos have left the system. Our simulations indicate that on 
average 90% of embryos in these systems were ejected during the course of the integration 
(i.e., either their semimajor axes exceeded 100 AU, or their orbital eccentricities became 
larger than unity) and among them, 60% collided with other protoplanetary bodies prior to 
their ejection from the system. A small fraction of embryos (~ 5%) also collided with the 
primary or secondary star, or with the Jupiter-like planet of the system. 



(1994 


), and 


David et al. 


(2003) 



thors, in binaries with small perihelia, the interactions of small bodies with the secondary star 
shorten the lifetimes of these objects and enhance the disk truncation. In binary-planetary 
syste ms, in addition to the pertu r bation from the stellar companion, similar to our solar sys- 



l2 



tem, (iChambers fc Cassen 1120021 : iLevison fc Agnor II2003I : [Raymond et al. II2004J : [Raymond 



20061 ). planetesimals and protoplanetary objects are also subject to the perturbative effect of 
the giant planet of the system. In such systems, the Jovian-type planet plays the important 
role of transferring angular momentum from the secondary star to planetary embryos and 
strongly affects the motion and radial mixing of these objects. Figure 6 shows this in more 
details. The systems simulated here are binaries with 0.5 solar-masses secondary stars and 
separations of 30 AU. The binary eccentricity in these systems is equal to 0, 0.2 and 0.4, from 
top to bottom. As shown here, as the eccentricity of the binary companion increases, its 
perihelion becomes smaller and its interaction with the giant planet becomes stronger. The 
latter causes the eccentricity of the giant body to increase and results in its closer approach 
to the disk of planetary embryos and enhancing collisions and mixing among these objects. 
The eccentricities of embryos, at distances close to the outer edge of the protoplanetary 
disk, rise to higher values until these bodies are ejected from the system. In binaries with 
smaller perihelia, the process of the transfer of angular momentum by means of the giant 
planet is stronger and the ejection of protoplanets occurs at earlier times. As a result, the 
total water to mass ratios of such systems become smaller as the eccentricities of their stellar 
companions increase. 

To further study the effect of the stellar companion on the dynamics and radial mixing 
of embryos, we carried out several simulations without the Jupiter-like planet of the system. 
Results indicate that, it is indeed possible to form terrestrial-class planets, with significant 
amounts of water, in the habitable zone of the primary star. However, because of the lack 
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of the intermediate effect of the Jovian-type planet, the interaction of embryos is slower 
and terrestrial planet formation takes longer. Figure 7 shows this for three systems. The 
separation of the binary in each system is 30 AU, and the mass of the secondary star is 1 
solar-mass. The eccentricity of the secondary star is equal to 0, 0.2, and 0.4 in simulations 
from left to right, respectively. 

An interesting result depicted by figure 7 is the decrease in the number of the final 
terrestrial planets, and increase in their sizes and accumulative water contents with increasing 
the eccentricity of the secondary star. The simulation on the left, in which the secondary is 
in a circular orbit, shows that, since in this system, the transfer of angular momentum from 
the stellar companion to protoplanetary objects, by means of the Jupiter-like planet of the 
system, is non-existent, the radial mixing of embryos is slow and inefficient. In binaries with 
larger eccentricities, the close approach of the stellar companion to the disk of protoplanets 
increases the rate of the interaction of these objects and enhances their collisions and radial 
mixing. As a result, in such systems, more of the water-carrying embryos participate in the 
formation of the final terrestrial planets. It is important to emphasize that, as explained 
below, this process is efficient only in moderately eccentric binaries. In binary systems with 



high eccentricities (small perihelia), embryos may be ejected from the system (jPavid et al. 



20031 ) . and terrestrial planet formation may become inefficient. 



An important result shown by figure 7 is the existence of a trend between the binary 
perihelion distance and the location of the outermost terrestrial planet. Figure 8 shows this 
for a set of different simulations. The top panel in this figure represents the semimajor axis of 
the outermost terrestrial planet, aout-, as a function of the binary eccentricity, e^. The bottom 
panel shows the ratio of this quantity to the perihelion distance of the binary, g?,. As shown 
here, simulations with no giant planet (shown in black) follow a clear trend: Terrestrial plan- 
ets only f orm interior to ro u ghly 0.19 times the binary perihelion distance. This has also been 



noted by iQuintana et al. I (120071 . see their figure 9) in their simulations of terrestrial planet 
formation in close binary star systems. The smallest binary perihelion that allows terrestrial 
planets to form outside the inner edge of the habitable z one (0.9 AU) in t h ese sy stems is sim- 



ply 0.9/0.19 = 4.7 AU, comparable to the estimate by IQuintana et al. I (120071 ). In binaries 
with no giant planets. Sun-like primaries with companions with perihelion distances smaller 
than approximately 5 AU are therefore not good candidates for habitable planet formation. 
It is, however, important to note that, because the stellar luminosity, and therefore the lo 



cation of the habitable zone , are sensitive to stellar mass (IKasting. Whitmire fc Reynolds 



19931 : [Raymond et al. 1120071 ). the minimum binary separation necessary to ensure habitable 



planet formation will vary significantly with the mass of the primary star. 

In simulations with giant planets, figure 8 indicates that terrestrial planets form closer- 
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in. The ratio aout/Oii in these systems varies between approximately 0.06 and 0.13, depending 
on the orbital separation of the two stars. The accretion process in such systems is more 
complicated since the giant planet's eccentricity and its ability to transfer angular momentum 
are largely regulated by the binary companion. 

4. Habitable Planet Formation 

Despite the stochasticity of the simulations and the large size of the parameter-space, 
many of our integrations resulted in the formation of Earth-sized objects, with substantial 
amounts of water, in the habitable zone of the primary star. Figures 9 and 10 show the 
results of some of these simulations. The orbital parameters of the final objects and their 
water contents arc listed in Table 1. It is important to mention that, in comparing the water 
contents of the Earth-like planets of our simulations with those of the Earth, we consider the 
water to mass ratio of the Earth to be 0.001. Since the exact amount of water in the Earth's 
mantle is unknown (between 1-10 Earth's ocean), such an estimate of Earth's water-mass 
fraction is equivalent to considering one ocean of water (~ 1.5 x 10^^ g) on the Earth's surface 
and three oceans in its mantle. 

A detailed analysis of the results depicted by figures 9 and 10 indicate that the systems 
shown in these figures have relatively large perihelia. Figure 11 shows this for simulations 
with = 1, in terms of the semimajor axis and eccentricity of the stellar companion. The 
circles in this figure represent those systems whose parameters were chosen from figure 1 (i.e., 
stable giant planets) and their simulations resulted in the formation of habitable bodies. The 
number associated with each circle corresponds to the mean eccentricity of the giant planet 
during the simulation. For the sake of comparison, the systems in which the giant planet is 
unstable have also been marked. Given that at the beginning of each simulation, the orbit 
of this planet was considered to be circular, a non-zero value for its average eccentricity is 
the result of its interaction with the secondary star. The fact that Earth-hke objects were 
formed in systems where the average eccentricity of the giant planet is small implies that 
this interaction has been weak. In other words, binaries with moderate to large perihelia 
and with giant planets on low eccentricity orbits are most favorable for habitable planet 
formation. Similar to the formation of habitable planets around single stars, where giant 
planets, in general, play destructive roles, a strong interaction between the secondary star 
and the giant planet in a binary-planetary system (i.e., a small binary perihelion) increases 
the orbital eccentricity of this object, and results in the removal of the terrestrial planet- 
forming materials from the system. 
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5. Summary and Conclusions 

We presented the results of the numerical simulations of the formation of Earth-like 
bodies in the habitable zones of moderately close binary-planetary system. The systems 
of our interest had binary separations equal to or smaller than 40 AU, their primary stars 
hosted Jupiter-like planets, and their habitable zones were closer to their primaries than 
the orbits of their Jupiter-like objects. We simulated the late stage of terrestrial planet 
formation in such systems, and studied the effect of the binary companion on the interaction 
of planetary embryos and the formation, dynamical evolution, and physical characteristics 
of the final terrestrial bodies. Our results indicate that, terrestrial-class planets can form in 
the habitable zones of binary-planetary systems, and their final assembly, as well as their 
water contents, strongly depend on the separation of the binary and the eccentricity of the 
stellar companion. 

Our simulations show that, similar to the formation of Earth-like objects at the presence 
of Jupiter in our solar system, where the ecce ntricity of the giant planet becomes the key 



factor in delivering water to the habitable zone (jChambers fc Cassen 1120021 : [Raymond et al. 



2004 ). the orbital dynamics of this object, and its intermediate role in transferring the 
perturbative effect of the secondary star to protoplanetary objects, in a binary planetary 
system, have strong effects on the radial mixing of embryos and the water-delivery process. 
As shown by [Raymond I (120071 ). giant planets are, in general, unfavorable for delivery of water 
from asteroidal regions to the habitable zone. In a binary system, the interaction between 
this object and the secondary star, particularly in eccentric binaries, increases its orbital 
eccentricity and makes the process of water delivery inefficient. Our calculations show that, 
to form habitable planets at the presence of a Jupiter-like planet in a binary star system, 
this interaction has to be weak. Wate r delivery is more efficient when the perihelion of the 
binary is large (|Quintana et al. 1120071 ) and the orbit of the giant planets is close to a circle 
jRavmond IboOQ l 



As mentioned in the introduction, this study was based on two fundamental assump- 
tions: (1) terrestrial planet formation in binary-planetary systems follows the same process 
as in single stars, (2) planetary embryos have already formed through oligarchic growth, 
and terrestrial-class objects are formed through the gravitational interaction of these bodies. 
Despite the observations of circumstellar disks around the components of several binary star 
systems, such assumptions may not be entirely viable. It is, in fact, necessary to study how 
the presence of a stellar companion affects the formation of planetary embryos, and also 
how the chemical structure of a disk of protoplanetary objects changes at the presence of a 
second close-by star. 



The final water contents of the terrestrial planets is also an issue that requires detailed 
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considerations. In our simulations, we assumed that all collisions were perfectly inelastic, 
and the water contents of the resulted planet would be equal to the sum of the water contents 
of the impacting bodies. This is an assumption that sets an upper limit for the water budget 
of terrestrial planets, and ignores the loss of water d ue to the impact. In fact, as shown by 
Genda fc Abe I (120051 ) and lCanup fc Pierazzo I (120061 ) . portions of the water of the impacting 
bodies will be lost due to the impact and the motion of the ground of an impacted body. 
The total water budget of the final bo dies of figures 9 and 1 may in fact be 5-10 times 
smaller than those reported in Table 1 (IRaymond et al. II2004I ). 



The simulations presented here are low resolution. Our model includes only large ob- 
jects. A more realistic model would contain smaller bodies such as km-sized planetesimals, 
as well. For this reason, the results of our integrations may not reveal detail characteristics 
of the fin al planetary systems, as they woul d be o btained from high resolutio n simulations 



[e.g., see [Raymond. Mandel 



fc Sigurdsson I (120061 )]. However, as shown by lAgnor et al. 



(119991 ) and [Chambers I (j200ll ) for the simulations of our solar system, such integrations can 
produce the main general properties of the final assembly of the planetary bodies. That, 
combined with the fact that the speed of computational simulations varies as N"^, where N 
is the number of involved objects, makes the low resolution simulations ideal for exploring 
the system's parameter-space. 

Despite of its limitations, our study shows that habitable planet-formation can be ef- 
ficient in moderately close binary-planetary systems. The favorable systems seem to be 
binaries with moderate to large perihelia and giant planets outside the habitable zones of 
their primary stars and in low eccentricity orbits. This points the fact that, to develop a 
more comprehensive understanding of the formation of habitable planets in such environ- 
ments, it is important to obtain a better knowledge of the formation of giant planets in these 
systems. 
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Fig. 1. — Graph of the critical semimajor axis of a planet in an equal- mass binary star 
system for three different values of the separation of the binary. 
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Fig. 2. — Radial distribution of original protoplanetary objects. 



- 18 - 



c 

CD 
O 
O 
LU 



O 
C 

o 
o 

LU 



0.4 

0.3 

0.2 

0.1 

0.0 
0.4 

0.3 

0.2 

0.1 

0.0 
0.4 

0.3 

0.2 

0.1 



L Myr J 


L 0.1 Myr J 

• • 

• 

. • ••••• 


L 1 Myr J 
• 

L • J 
= ^ '• *• • • • • 

: 1 1 1 i-^-^-^-^-n^^^K-,; 


^^^W . 

L 10 Myr J 

. • 


L ^ ^ 30 Myr J 

L • J 
• 


: 1 1 1 1 1 i i ■: 

L 100 Myr J 

:•• 





1 


2 3 4 5 1 
Semimajor Axis (AU) 

Log(Water Mass Fraction) 


2 3 4 
Semimajor Axis (AU) 


5 






1 


\ ^^^^^H 




1x10 ° 




1 






-5 




-4 -3 


-2 


-1.3 



Fig. 3. — Snap shots of the interaction of protoplanetary objects and the formation of 
terrestrial planets. The mass of the secondary star is 0.5 solar-masses and its semimajor axis 
and eccentricity are 30 AU and 0.2, respectively. The results show a terrestrial-sized planet, 
with substantial amounts of water, at a semimajor axis of 1.2 AU and with an eccentricity 
of approximately 0.07. The Jupiter-sized planet of the system is shown by the black circle. 



- 19 - 



0.4 



i 0.2 

o 

o 

LU 

0.1 



0.0 
0.4 



i 0.2 

o 

o 

LU 

0.1 



0.0 
0.4 



S 0.2 

o 

o 

LU 

0.1 



0.0 



L 1 00 Myr j 
• 

fi 1 1 1 1 


L 100 Myr J 

.. • J 


L 1 00 Myr j 
• 

; • 


L 100 Myr J 
• 


L 1 00 Myr j 
• 


L 100 Myr J 

r • 



2 3 4 

Semimajor Axis (AU) 



12 3 4 

Semimajor Axis (AU) 



Log(Water Mass Fraction) 



Fig. 4. — Stochasticity of simulations. The left column depicts the results of three simulations 
with random distribution of protoplanetary objects in a binary system with a Sun-like star 
as its secondary and in a circular orbit with a radius of 30 AU. The right column shows the 
results of simulations for similar distributions of planetary embryos in a system with a 1.0 
solar-mass secondary star in an orbit with a semimajor axis of 40 AU and eccentricity of 
0.2. As shown in each column, different distribution of planetary embryos may result in the 
formation of different number of planets with different water to mass ratios. 



- 20 - 



a =40 AU 
' e =0.4 

b 



a =40 AU 
' e =0.2 

b 



a =40 AU 

' e=0 

b 



a =30 AU 
' e =0.2 

b 



a =30 AU 

' e=0 

b 



a =20 AU 

' e=0 

b 



Solar 
System 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



HZ 



I I I 



I I I 



• • • 

I I I 



• 9 



• •- 



• • • 

LU LJJ I I I 



Lil LU LU 



• • • 

LU iij m 



I I I 



I I I 



LU 



I I I 



I I I 



I I I 



I I I I I I I I I 

1 2 3 4 5 

Semimajor Axis (AU) 

Log(Water Mass Fraction) 



5X10 
1x10° 



Fig. 5. — Results of simulations in a system with /if, = 1 for different values of the eccentricity 
(cfo) and semimajor axis (at) of the stellar companion. 
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Fig. 6. — Variation of water contents with the eccentricity of the stellar companion. The 
primary star has a mass of 0.5 solar-masses, the semimajor axis of the binary is 30 AU, and 
its eccentricity is equal to 0,0.2, and 0.4, from top to bottom. 
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Fig. 7. — Results of simulations in binary systems with no Jupiter-like planet. The stars 
of each binary are Sun-like and their separations are 30 AU. The orbital eccentricity of the 
secondary star is 0, 0.2, and 0.4, for the systems on the left, middle, and right, respectively. 
Note that the time of integration has been increase to 200 Myr. 
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Fig. 8. — Top panel: Semimajor axis of the outermost terrestrial planet. Bottom panel: 
The ratio of the semimajor axis of this object to the perihelion distance of the binary. The 
secondary star is solar mass. 
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Fig. 9. — Earth-like objects in the habitable zone of the primary star. The mass of the 
secondary star in all simulations is 1 solar-mass. 
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Fig. 10. — Habitable planet formation in binary-planetary systems with 0.5 and 1.5 solar- 
masses secondary stars. 
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Fig. 11. — Habitable planet formation in the (e;,, a^) space of an equal- mass binary-planetary 
system. Circles correspond to binaries with initial parameters chosen from figure 1, in 
which habitable planets are formed. Triangles represent systems in which the giant planet 
is unstable. The number associated with each circle represents the mean eccentricity of the 
giant planet of the system during the simulation. 
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Table 1. Characteristics of the Earth-Uke planets of figures 9 and 10. 
(Sun's habitable zone: 0.9 - 1.50 AU) 



Simulation 


Planet Mass (M®) 


Semimajor Axis (AU) 


Eccentricity 


Water/Mass 


8-A 


0.95 


1.28 


0.03 


0.00421 


8-B 


0.75 


1.11 


0.06 


0.00415 


8-C 


1.17 


1.16 


0.03 


0.00164 


8-D 


0.86 


1.33 


0.09 


0.01070 


8-E 


0.95 


1.50 


0.08 


0.00868 


9-A 


0.74 


1.07 


0.06 


0.00349 


9-B 


0.99 


1.26 


0.12 


0.00366 


9-C 


1.23 


1.30 


0.09 


0.00103 



